The Annals of Probability 
2008, Vol. 36, No. 2, 765-795 
DOI: 10.1214/07-AOP343 

© Institute of Mathematical Statistics, 2008 



ON THE ORTHOGONAL POLYNOMIALS ASSOCIATED 
WITH A LEVY PROCESS 1 

By Josep Llui's Sole and Frederic Utzet 



Let X — {X t , t > 0} be a cadlag Levy process, centered, with mo- 
ments of all orders. There are two families of orthogonal polynomials 
associated with X. On one hand, the Kailath-Segall formula gives the 
relationship between the iterated integrals and the variations of order 
n of X, and defines a family of polynomials Pi(a;i), P2(xi, X2), ■ ■ ■ that 
are orthogonal with respect to the joint law of the variations of X. On 
the other hand, we can construct a sequence of orthogonal polyno- 
mials Pn(x) with respect to the measure a 2 6o(dx) + x 2 v(dx), where 
a 2 is the variance of the Gaussian part of X and v its Levy measure. 
These polynomials are the building blocks of a kind of chaotic rep- 
resentation of the square functionals of the Levy process proved by 
Nualart and Schoutens. The main objective of this work is to study 
the probabilistic properties and the relationship of the two families 
of polynomials. In particular, the Levy processes such that the asso- 
ciated polynomials P n (xi, . . . ,x n ) depend on a fixed number of vari- 
ables are characterized. Also, we give a sequence of Levy processes 
that converge in the Skorohod topology to X, such that all variations 
and iterated integrals of the sequence converge to the variations and 
iterated integrals of X. 



1. Introduction. Let X = {Xf,t > 0} be a semimartingale with Xq = 0. 
Define the iterated integrals by the recurrence 
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and consider the sequence of the variations of X, 

(1.2) Af ) = X t , X[ 2) = [X,X] t , X™ = ( AX sT, n>3, 

0<s<t 

where AX S = X s — X s _. The Kailath-Segall formula (see Segall and Kailath 
[7] or Meyer [12]) gives the relationship between P^ and x\ n ^ : 

(1.3) p t {n) = -(p^x^ - p^x[ 2) + ■■■+ (-ir +i p, (o) x, (n) ). 

n 

We deduce that P^ is a polynomial in Xjp , . . . , xf 1 ' , called the Kailath- 
Segall polynomial of order n. Denote this polynomial by P n {x\, . . . ,x n ), so 

p^ = p n (x^,...,xl n) ). 

The explicit expression of P n (xi, . . . , x n ) is 

p n (x 1 ,...,x n )=(-ir5:n L ^ I T 5 

where the summation is over all nonnegative integers mi, . . . ,m n such that 
Y^j=i3 m 3 = n ( see Avram and Taqqu [3]). The polynomials P n {x\, . . . ,x n ), 
n > 1, are also a particular case of generalized Appell polynomials (see 
Anshelevich [1]). The first three of these polynomials are 

Pl{xx) = XI, 

P 2 (x 1 ,X2) = \x\ - \x 2 , 
P?,(xi,X2iXz) = ~x\ - \x\X-i + ~Xs. 

If X is a martingale with predictable quadratic variation (X,X) t = Ct, 
with finite moments of all orders, then p( n ) and X^ also have moments of 
all orders and the iterated integrals of different order are orthogonal, that 
is, 



(i.4) npt ) pl m) ]=-,c n t n 5 nn , 

n\ 

where 5 nm = 1 if n = m and otherwise, and C = E[Xf]. The orthogonality 

(2) 

of the Kailath-Segall polynomials with respect to the law of (Xt,X\ , . . .) 
follows. This is true, in particular, for a centered Levy process with moments 
of all orders. 

Consider a centered Levy process X with moments of all orders and let a 2 
be the variance of its Gaussian part and v its Levy measure. The measure 
^"{dx) = a 2 5o(dx) + x 2 v(dx) [we also write ^{dx) = x 2 is(dx)] is finite and 
we can construct a (finite or infinite) sequence of orthogonal polynomials 
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Pni x ) [resp. p n (x)] with respect to 7 " (resp. 7). These determine a sequence 
of strongly orthogonal martingales related to the Teugels martingales (see 
Nualart and Schoutens [13]) that are the building blocks of a kind of chaotic 
representation of the square functionals of the Levy process. Therefore, we 
will call {p°(x),n > 0} [or p n (x), when a = 0] the Teugels polynomials as- 
sociated with X. 

The main objective of this paper is to study the probabilistic properties 
of, and the relationship between, the polynomials P n (xi, ■ . ■ , x n ) and Pn(x). 
When working on that problem, we found three results that we think are 
interesting in themselves. The first one is a proof that for a general semi- 
martingale X, the Doleans exponential £{uXt) (fixed t and uj) is analytic 
in a certain neighborhood of the origin and that the iterated integrals are 
the Taylor coefficients. This part is based on a paper of Lin [11], where the 
result for a Levy process was implicitly proven, and on Yablonski [19] , where 
a generating function of the Kailath-Segall polynomials of a Levy process 
was introduced. 

The second result is related to the Kailath-Segall polynomials that are 
expressible as polynomials of a fixed set of variables. A very interesting 
property of the Kailath-Segall polynomials is that when you impose some 

restriction on the variables xi,X2,.-., you get different well-known families 

(2) 

of polynomials. For example, when A" is a Brownian motion, then X f = t 

(n) 

and XI = 0,n > 3, showing that it is enough to consider the polynomials 
P n (x, t, 0, . . . , 0), and it turns out that 

P n (x,t,0,---,0) = H n {x,t), 

where H n (x,t) are the (generalized) Hermite polynomials. In a similar way, 
considering a compensated Poisson process, the (generalized) Charlier poly- 
nomials C n (x,t) are obtained. It is known that the Brownian motion and 
the compensated Poisson process are the unique Levy processes such that 
the Kailath-Segall polynomials (that means, the iterated integrals) can be 
written as polynomials in x and t (see Section 3). So, a natural question is 
how to characterize the Levy processes with a similar property for a finite 
number of variables. The answer is that they are the Levy processes such 
that the Levy measure has finite support. This is not surprising, given the 
paper of Sengupta and Darkar [16], where a similar result was obtained in 
relation to space-time harmonic polynomials. The key to our proof is that, 
under the appropriate conditions, only the application of linear functions to 
a Levy process gives rise to another Levy process 

Finally, the third result that we would like to mention is that, under 
the appropriate hypothesis, it is possible to give a sequence of simple Levy 
processes {X^,k > 1} that converges in the Skorohod topology to X and 
these processes satisfy the conditions of Avram [2] in order that all variations 
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and iterated integrals of converge to the variations and iterated integral 
of X. This approximating sequence is constructed using the Gauss-Jacobi 
mechanical quadrature formula. 

2. Doleans exponential and Kailath Segall polynomials. This section is 
inspired by the works of Lin [11] and Yablonski [19]. First, in the paper of 
Lin [11], it is implicit that the iterated integrals of a Levy process are the 
Taylor coefficients of the Doleans exponential at the origin, a property sug- 
gested by Meyer [12], page 318, for a semimartingale. Second, Yablonski [19] 
introduced a generating function in order to study a family of polynomials 
associated with a Levy process that turn out to be the Kailath-Segall poly- 
nomials. However, that generating function is deterministic and Yablonski 
gives no probabilistic interpretation of it. Here, we combine both approaches 
for a general semimartingale and prove that the Yablonski generating func- 
tion is the Doleans exponential of the semimartingale for fixed u>, which is 
analytic in a neighborhood of the origin, and that the Taylor coefficients are 
the iterated integrals. Therefore, we prove the general claim of Meyer [12]. 

To begin with, for the sake of easy reference, in the next remark, we collect 
some results obtained by Yablonski [19]. 

Remark 2.1. Given a sequence of real numbers x = (x\,X2, ■ ■ ■) such 
that limsup n l^nl 1 ^™ = A < oo, Yablonski [19] defines the generating function 

r oo /-_ 1 \n+i 
F(u,x) = exp< ^ u n x, 

ln=l ^ 

which is analytic for u G ( — 1/A, 1/A). Yablonski proves that in the expansion 

oo 

(2.1) F(u,x) = J2u n Pn(x), 

n=0 

the function P n (x) is a polynomial in x±, . . . ,x n , which satisfies 

Pn {%1 j • ■ • j ^-n) 

= -(P n -i(xi, ■ ■ .,x n -i)xi 
n 

- P n - 2 (xi,. . .,X n - 2 )x 2 H h (-l) n+1 Po^n), 

where Po = 1. Comparing this with (1.3), we deduce that P n = P n . 

Further, Yablonski [19] also points out the following very useful properties: 

P n (flXi , . . . , (L Xn) — CL P n (xi , • • • , Xn) 

and 

n 

(2.2) P n (xi +yx, ...,x n + y n ) = Y^ Pk(xi,-- ■ , x k )P n -k(yi, ■ ■ -lUn-h)- 

k=0 
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Let X = {Xt,t > 0} be a semimartingale with Xq = 0. For uGl, consider 
the Doleans equation 

Z t = l + u f Z s -dX s , 
Jo 

which has a unique solution (semimartingale) given by the Doleans expo- 
nential of uXt, 

(2.3) S{uX t ) = exp{uX t -±u 2 (X c ,X c ) t } [J (1 + uAX s )e~ uAXs . 

0<s<t 

Fixing uj and t € R+, it is clear that if {X s ,s € [0,t]} is continuous or has a 
finite number of jumps, then £{uXt) is analytic for u£l. The proposition 
below provides a general result in this direction. 



Proposition 2.2. Fix wef! (out of a set of probability zero) and t G 
There then exists uq>1, depending on ui and t, such that the function 
£{uXt) is analytic in u E (— uq,uq) and 



(2.4) £{uX t ) = Y,u n Pr- 



oo 

n n(") 



n=0 



Proof. Fix us € 0. From the expression (2.3), it follows that we only 
need to prove that Ilo<s<t(l + uAX s )e~ uAXs is analytic. Decompose this 
product in the following way: 



[ (l + uAX s )e- ul 



>j,AX 3 
0<s<t 

-uAX. 



J] (l + uAX s )e- uAX ° [] (l+«AI s )e- 

0<s<t 0<s<t 
|AX S |<1 |AJf s |>l 



(*) (**) 

The term (**) is analytic, since there is only a finite number of factors. On 
the other hand, for u G (—1,1), the expression (*) is positive and, taking 
logarithms, we have 

0<s<t n=2 
|AX S |<1 

= wE LJ ^- n E (A*.r, 

n=2 0<s<t 
|AXJ<1 
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where (a) follows from log(l + y) — y = Y^=2 — 2/ n > the series being 



absolutely convergent for \y\ < 1 and (b) is due to Fubini's theorem, which 
can be applied since 

n 



E E 

0<s<t n=2 
\AX S \<1 



n 



-u n (AX s ) n 



00 \n\ n 

E^- E lAx.r 

n=2 0<s<t 
|AX S |<1 



<C (&X S ) 2 <C[X,X] t <oo. 

0<s<t 
|AX S |<1 



In a similar way, it is computed that 

l/n 



lim sup 



\n+l 



E ( A ^) n 

0<s<< 
|AX S |<1 



/ 1 \ 1 / n 

<limsupf -[X,X] t J <1. 



The analyticity of (*) for u G ( — 1,1) then follows. 

To compute the coefficients of the expansion of £(uXt), let u\ = 
(maxo<s<t |AX S |) _1 . Then, for u G (— |nAX s j < 1, Vs G (0, t]. There- 
fore, we can repeat the preceding proof to obtain 

Y[ {l+uAX s )e- uAX ° 

0<s<t 



(2.5) 



exp -± £ - 2 (AX S ) 2 + £ £ 



0<s<i 



0<s<tn=3 



and 



lim sup 



.-nn+i 



E ( A ^) r 

0<s<< 



l/n 



< lim sup 



U G 



l/n J 

< 

«1 



Since the right-hand side of (2.5) is the generating function of Yablonski, 
from Remark 2.1 we deduce the expression (2.4) for u G (— u\,u\). If u\ < 1, 
by the principle of analytic continuation, we deduce (2.4) for u G (—1,1). 
Finally, we take uq = max{l, u±}. □ 



Remark 2.3. From the preceding proof, we deduce that for u£ (—1, 1), 



£(uXt) = exp(nX 4 - \u 2 {X c ,X c ) t 



+ / (log(l + ux) — ux) dJ(s, x) 

'(0,*]x{0<|x|<l} 
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x n (l + uAX s )e- uAX % 

0<s<t 
|AX S |>1 

where 

J(B) = #{t : (t, AX t ) eB}, Be B((0, oo) x M ), 

is the jump measure of X, where Mo = ^ ~~ {0} (see Jacod and Shiryaev 
[6], page 69). Moreover, fixing uj G VL and t > 0, for -u € (— iti,Ui) where 
u\ = (maxo< s <t (AA^I)" 1 (which depends on u and t), 

£{uX t ) = expj p t^ u n x (n) I 

= exp\uX t --u 2 (X c ,X c ) t + [ (log(l + ux) - ux) dJ (s , x)\ . 
I 2 J(0,t]xR J 

For any semimartingale X, we have that (aX)^ = a n X( n ). Further, given 
two semimartingales X and Y such that [X, Y] = 0, we have 

(I + 7)( n )=l('')+y( n ) Vn>l. 

This can be proven from the expression for the product of two Doleans 
exponentials, 

S(uX t )£(uY t ) = 8{uX t + uY t + u 2 [X, Y] t ), 

formula (2.2) and formula (2.4). We summarize these formulas in the next 
proposition. 

Proposition 2.4. Let X and Y be two semimartingales such that Xq = 
Y = and [X, Y] = 0. For a,beR, write Z = aX + bY . Then, 

P n (Z?\...,Z^) 

= £ aH n - k P k (X?\. . . , Af y B _ fc (y t (1) , . . . , Y t ^ k) ). 

k=0 

3. Kailath Segall polynomials associated with a Levy process. From now 
on, consider that X = {Xt,t > 0} is a Levy process (meaning that X has 
stationary and independent increments, is continuous in probability and has 
Xq = 0), cadlag, centered and with E[|Xi| n ] < oo for every n > 1. Denote by 
a 2 the variance of the Gaussian part of X and by v its Levy measure. Since 
it is a martingale with predictable quadratic variation 



(X,X) t = {a 2 + j^x 2 u{dx)^t, 



8 J. L. SOLE AND F. UTZET 

it follows that the Kailath-Segall polynomials are orthogonal. 
When X is a Brownian motion W = {Wt,t > 0}, 

=W t , X[ 2) =t and X (n) =0, n > 3, 

and, therefore, 

P n (X[ l) , . . . , X[ n) ) = P n iW t , t , 0, . . . , 0) = Q n (W t , t) . 

That is, it is enough to consider the Kailath-Segall polynomials with the 
variables 

x\ = x, X2 = t and x n = 0, n > 3, 

and then 

P n (x,t,0, ...,0) = Q n (x,t) =H n (x,t), 

where H n {x,t) are the Hermite polynomials defined via the generating func- 
tion 

oo 

exp-fux — \u 2 t} = ^ u n H n (x,t). 

n=0 

Note that the leading coefficient of H n (x,t) is 1/n!, and for fixed t > 0, 
H n (x,t) and H m (x,t), n^m, are orthogonal with respect to the Gaussian 
measure M(0,t). However, observe that H n (x,t) is an ordinary polynomial 
in x and t and is defined for all For t = 0, 

P n (x,0,0,...)=#nOr,0) = ^. 

n! 

If X is a compensated Poisson process of parameter b > and jumps size 
a, that is, Xt = a(Nt — bt), where N = {Nt,t > 0} is a standard Poisson 
process of intensity b, then 

X[ n) = a n N t , n>2. 

Therefore, 

P n {x[ l) ,...,X?) = P n (X t ,aX t + a 2 bt, a n ~ x X t + a n bt) = Q n (X t ,t) 

and the polynomial Q n (x, t) can be explicitly computed in the following way: 
write x = a(y — bt), then 

Q n (a(y-bt),t) = a n C n (y,bt), 

where C n (x,t) is the Charlier polynomial with leading coefficient 1/n! de- 
fined by 

oo 

e- tu (l + u) x = Y,u n C n (x,t). 

n=0 
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Again, note that C n (x,t) is defined for every x,t € M and, in particular, 

C n (x , 0) = P n (x , x , . . . ) = — j— , 

n! 

where [x] n is the falling factorial, [x] n = x(x — 1) • ■ • (x — n + 1), [rr]o = 1. 
Fixing £ > 0, the polynomials C n (x,t),n > 1, are orthogonal with respect to 
the Poisson distribution of parameter t. 

Moreover, it is known that the Brownian motion and the compensated 
Poisson process are the unique Levy processes such that the Kailath-Segall 
polynomials can be written as polynomials in x and t. This fact follows from 
Feinsilver [8], who gives a necessary condition for the iterated integral 
to be a polynomial in Xt, that condition being satisfied only by the binomial, 
negative binomial, Gamma, Poisson and Gaussian types (see Feinsilver [8], 

(2) 

page 301). It is easy to check that P t is not a polynomial on Xt for the 
binomial, negative binomial and Gamma process; see also Privault et al. [14] 
for a different proof. 

Therefore, we can ask if there are Levy processes such that the Kailath- 
Segall polynomials depend on a fixed finite set of variables. The answer is 
affirmative and the examples are very easy to find. From Proposition 2.4, 
we deduce, for the process 

X t = aWt + a(N t -bt), 

where W is a Brownian motion and N is a Poisson process of parameter 
b > 0, independent of W, that the Kailath-Segall polynomials can be written 
as polynomials in yo, y\ and t, which are the convolutions of the polynomials 
aH.(yo,t) and aC.{y\,bt) described above. 

More generally, a jump diffusion Levy process 

X t = aW t + J t , 

where J is a centered compound Poisson process with only a finite number 
of jump sizes, has Kailath-Segall polynomials expressible in a fixed finite 
set of variables. Specifically, let 

n 

X t = aWt + Y,a j {N j {t)-b j t\ 

3=1 

where W is a Brownian motion, Nj is a Poisson process of parameter bj, the 
processes W, N\, . . . , N n are independent and a±, . . . , a n are different nonzero 
numbers. These kinds of processes will be called simple Levy processes and 
will play a key role. For such processes, we will see that there is a family of 
polynomials Q m (xi, . . . , x n +2) such that, for m > n + 2, 
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Moreover, for m > n + 2, let 7Z™ be the subspace of M. m given by the vectors 
(xi, . . . , x m ) such that there exists (yo, . . . , y n ) € M n+1 with 

n 

x\ = vyo + ^ajivj - bjt), 

3=1 
n 

x 2 = a 2 t + ^2a 2 j y j , 

3=1 

n 

x k = ^a k j y j , k = 3,...,m. 
i=i 

Then, the polynomial P m restricted to H™ is a (multiple) convolution of 
aH.(y ,t) and a,jC.(yj,bjt), j = l,...,n. 
If a = 0, then 

P m {X[ l) ,...,xt ] ) = Qm (X[ l) ,...,X^ n \t), m > n + 1, 

and in the expression for Q m restricted to a similar subspace as above, there 
is no Hermite polynomial part. 

Remark 3.1. To summarize the situation, P n (x±, . . . ,x n ),n > 1, is a 
family of ordinary polynomials that can be evaluated on an arbitrary se- 
quence of real numbers or random variables. However, the most interesting 
properties of P n (x\, . . . ,x n ) appear when we consider a centered Levy pro- 
cess with finite moments of all orders and apply P n on the sequence of the 
variations of X : (X^\ X^ 2 \ . . .). Then, by the Kailath-Segall formula (1.3), 

P^prf,..,!*" 1 ), 

and P n {X^\ . . .,x[ n) ) and P m {x[ l \. . . ,X t (m) ) are orthogonal if n± m. 

From an equivalent point of view, let M°° = {(x±,X2, ■ ■ -),x n € M} be the 
set of sequences of real numbers and the probability on (R°°, ^(IR) 00 ) 

induced by {x[ l \xf\. . .). Writing P n {xi, ...,x m ,...) = P n (xi, . . . ,x n ), P n 
can be considered as a polynomial defined on M°° and the different polyno- 
mials are orthogonal: 



/ P n P m dP t °° = if n ^ m. 



In some cases, the probability is concentrated in a finite-dimensional 
subspace of M°° and then the restriction of P n to this subspace gives rise to 
a new family of polynomials that depend on a finite set of variables. 

So, a natural question is whether there are other examples, different from 
simple Levy processes, where the Kailath-Segall polynomials can be written 
as polynomials in a finite number of variables. We will prove that the answer 
is "no." 
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3.1. Polynomials of a Levy process. The purpose of this subsection is 
to study when a polynomial of a Levy process can be a Levy process. The 
result is given in the following proposition. 

Proposition 3.2. Let {(Y(i),Xi(i), . . . ,X d (t)),t > 0} be a d+l- dimen- 
sional Levy process with moments of all orders and P(x\, . . . ,Xd,t) a poly- 
nomial. IfY(t) = P(Xi(t),...,Xd(t),t), then P(x\, ... ,Xd,t) has degree 1. 

In order to prove this proposition, we will need the following elementary 
property. 

Lemma 3.3. Let A be a dx d nonnegative definite matrix o/rang(A) = 
r <d and let f = (/i(t), • ■ • , fd{t))' be a vector of real functions such that 

f'Af = 0. 

Then: 

1. if r = d, then f = 0; 

2. if r < d, then there are d — r functions between fi,-.-,fd, such that the 
other r depend linearly on them. 

Proof. This is a consequence of the fact that for any matrix C con- 
formable with A such that CAC = 0, we have AC = 0. □ 

PROOF of Proposition 3.2. The idea of the proof is that from the Ito 
formula, the decomposition of If as a special semimartingale is obtained, 
and Jacod and Shiryayev [6], Corollary II. 4. 19, give necessary and sufficient 
conditions in order that a semimartingale be a Levy process. Then, we will 
prove that for a polynomial of degree n of a Levy process to be a Levy 
process, it is necessary that a polynomial of order n — 1 be a Levy process. 
Hence, we can reduce to the case where the polynomial has degree 2. The 
proof is as follows. 

First, since a polynomial in (Xi(t), . . . ,Xd(t),t) can be written as a poly- 
nomial in (Xi(t) — K[Xi(t)], . . . ,Xd(t) — K[Xd(t)], t), we can assume that the 
Levy process is centered. Also, every linear combination Y^ = x\jXj{t) + fit 
is a Levy process (jointly with Y) and we can then eliminate such linear 
combinations from P(X\(t), . . . ,Xd(t),t); that is, we will assume that every 
monomial in P(x\, . . . ,Xd,t) has degree > 2. 

1. Degree of P(x±, . . . , x d , t) = 2. In this step, we will prove that if a 
quadratic form P(Xi(t), . . . , Xd(t),t) is a Levy process, then P(x±, . . . , x d , t) 
I). For now, we write Xd+i = t. Then, 

d+l 

P(xi , . . . , X d +l) = hiX< i + Yl 

CijXiXj . 

i=l i<j 
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First, by a linear transformation, we can assume that 

d+l 

P(x 1 ,...,x d+1 ) = ^2bixf. 



i=l 



We then have 



Y t = Y,bjX?(t) + b d+1 t 2 . 

Taking expectations, and recalling that the Levy process is centered, we get 

Ct = b d+1 t 2 . 

So, bd+i = 0. We will prove by induction over d that also b\ = ■ ■ ■ = b d = 0. 

1.1. Let d = 1. In order to prove that X 2 cannot be a Levy process, write 
Y t = X 2 and let C2 = E [X 2 ] = a 2 + J K x 2 v{dx) be the variance and cumulant 
of order 2 of X\ and C4 = j- R x i v{dx) be the cumulant of order 4 of X\. On 
one side, if Y were a Levy process, then E[(Yf — EYf) 2 ] = Ct, for some C . 
On the other side, from the relationship between moments and cumulants, 
we have 

E[Y t 2 } = E[Xf] = C±t + 3C 2 t 2 . 

Comparing both expressions for E[Y^ 2 ], we deduce that X = 0. 

1.2. Consider d > 2 and let X t = (Xi(i), . . . ,X d (t))' be given by 

X t = B 4 + / xN(ds,dx), 

where B 4 = (Bi(t), . . . ,B d (t))' is a (i-dimensional Brownian motion with co- 
variance matrix A and N(i,x) is the jump measure of the process, where 
Mq = R rf — {0} and <iN(i,x) = dN(t,x) — dtdv{-x) is the compensated jump 
measure. 

By Ito's formula, 

P(X t ) = 2 /VMQOHdB^) 

, Jo i=i 

(3.1) 

P d 

+ / V(26iXiXi(s-) + bixf) dN(ds, dx) 
J(o,t]x^ i=1 

(3.2) + 1 (j^ A abi + J Rd ^( dx )) • 

The right-hand side of (3.1) is a martingale and (3.2) is of bounded variation 
and continuous, so the above expression is the decomposition of P(X-t) as a 
special semimartingale. 
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1.2.1. Assume By Jacod and Shiryavev [6], Corollary II. 4. 19, a 

necessary condition for P (X-t) to be a Levy process is that the quadratic 
variation of the continuous martingale part should be of the form Ct (the 
truncation function does not play any role in that condition). Then, 



This implies that 



J2 kbjAijXi (t)Xj (t) = C Vt > 0, 



and from the fact Xj(0) = for i = 1, . . . , d, we deduce that C = 0. We can 
write this expression in a vector form, 

U' t AU t = 0, 

where XJt = (b\Xi(t), . . . , bdXdit))' . By Lemma 3.3, we see that if some bi ^ 
0, then there is a linear relationship between X±, . . . ,Xd- 
1.2.2. If A = 0, then 

„ d 

/ Y(2b i x i X i (s-) + b 2 i xf)dN(ds,d^) 

J(o,t]xB$ i=1 

is a Levy process. From the fact that J^ t ] xK d x\ dN(ds, fix) = [Xi,Xi]t, and 
as we are assuming that (Yt,Xt) is a Levy process, we deduce that 

„ d 

is a Levy process and a martingale. Therefore, 

V r d x 21 

E 



y^(biXiXi(s))dN(ds,dx) 
(0,t]xR^ 



E 



(0,t]xEg 



dsv(dx), 



J2(biXiXi(s)) 

\i=l / J 

but the left-hand side has the form Ct, and the right-hand side C't 2 . Hence, 
d 

*^2biXiXi(t) = y(x,t,oj), v{dx) (g> dt <g> P a.e. 



i=l 



From the fact that v ^ 0, it follows that there are x\, . . . ,Xd, not all 0, such 
that 

d 

^biXiXiit) = V(i,w), dt <g>P a.e. 
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and this implies a linear relationship between X\, . . . ,Xd- 

1.2.3. From 1.2.1 and 1.2.2, we then deduce that if P ^ 0, there is a 

linear relationship between X\, . . . and it follows that there is a Levy 

process (Y(t),X±(t), . . . and a polynomial of degree 2 such that 

Yt = P(Xi(t), . . . ,Xd,—i(t)). Iterating the procedure, we arrive at the case 

d = l, which is absurd (Point 1.1). 

2. Degree of P(x\, . . . ,Xd,t) = n > 3. This proof is very similar to 1.2.1 

and 1.2.2. As in point 1.2, we apply Ito's formula. 

2.1. Assume that the covariance of the Gaussian part of X is not zero: 

A ^ 0. The continuous martingale in -P(X^) is. 

1 dP 



rt " y>u 



By Jacod-Shiryaev [6], its quadratic variation should be Ct. Write V = 
(|£(X S , s), . . . , §f d (X s ,s))'. We then have 

V'AV = 0. 

By Lemma 3.3, it follows that there are numbers gi,... ,gd, not all null, such 
that 

d Qp 

i=l % 

and the expression on the left-hand side is a polynomial of degree n — 1 on 
Xx(t),...,X d (t),t. 

2.2. If A = 0, then consider the bounded variation part of P(Xf,t), 



^^'(^ (x *• s)+ ^( p(x *- +x ' s, 



d dP \ \ 

-P(X s _,s)-^:ri — (X s -,s)jv(dx)J ds, 

which is continuous, thus predictable. On the other hand, the Levy-Ito ex- 
pression for a Levy process also gives its decomposition as a special semi- 
martingale, so, by the unicity of the decomposition, we deduce that Vt = Ct. 
Then, 



dP 
~dt 



(Xt,t)+ / (p(X t -+x,t)-P(Xt-,t)-Y / x i ^(X t -,t))v(dx)=C 

This expression is also a polynomial of degree n — 1 in X\(t), . . . ,Xd(t),t. 

2.3. From 2.1 and 2.2, we deduce that for a polynomial of order n in 
X\(t), . . . ,Xd(t),t, to be a Levy process, it is necessary that a polynomial of 
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order n — 1 be a Levy process. Iterating, we arrive at a contradiction with 
step 2. □ 

Remark 3.4. An indication that the property expressed in Proposition 
3.2 may be true for more general functions is the following. Instead of a 
polynomial, consider a general (sufficiently regular) function f(x\, . . . ,x,j). 
Assume that the covariance matrix A of the Gaussian part of X is nonsin- 
gular. The necessary condition of Jacod and Shiryavev [6], Corollary II. 4. 19, 
for /(X t ) to be a Levy process becomes 

Xj(^( X *-)) =C a.s., for alii > 0, 

where g is a function obtained from / through linear changes of variable. 
Since the support of Xj is M. d , it follows that 

\\Vg(x 1 ,...,x d )\\ 2 = C y(x 1 ,...,x d )eR d . 

This is the eikonal equation, which, in M. d , has a unique solution given by a 
linear function; see Khavinson [9], Remark (ii), or Letac and Pradines [10]. 

3.2. The n -dimensional variation process We return to 

the general Levy process X with moments of all orders. From E[|Ai| fc ] <oo, 
for all k > 1, it follows that E[|A^ n) | fc ] < oo, for all n,k>l and 

E[A f (1) ] = 0, E[A t (2) ] = (a 2 + J^ x 2 u{dx)^j t 

and 

E[X t {n) ] =t f x n v(dx), n>3. 
Jr 

Consider the multivariate Levy process (X^ , . . . , x[ n ^). Its Levy measure 
u n (on W 1 ) is the image measure of v by the application 

E — >R n , 

x i— > (x, x 2 , . . . , x n ). 

By the image measure theorem, for / : M n — ► R measurable, positive or u n - 
integrable, 

/ /(x)i/ n (dx) = / f{x,x 2 ,...,x n )v(dx). 

JM. n Jr 

The characteristic function of (x\ l \ . . . ,X^ ) is then 

(3.3) 99 t (z) = exp|-itzfcj 2 + itz 2 a 2 +t J (e i ^i= 1 * iX ' -l-izix)u(da)\, 
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where z = (z±, . . . , z n ). Hence, the characteristic function of a linear combi- 
nation 

n 

Zt = ^ a jXt^ + On+lt, 
3=1 

such that E[Z t ] = 0, is 

4>t{z) =exp< — ^ta\z 2 a 2 

(3.4) 



v{dx) > , z 6 ' 



Before proceeding to the main theorem of this section, we need the fol- 
lowing lemma which will allow us to work with characteristic functions like 
(3.4). 

Lemma 3.5. Let v be a Levy measure on (R,23(R)) and let /:R— > M 6e 
a continuous function such that /(0) = and f^f 2 {x)v{dx) < oo. // 

( e te /(*)-l-fe/(x))i/(dx) = VzGK, 



then f = 0, v-a.e. 

Proof. Let z/j be the measure image of v by /. From the hypothesis, 
it follows that Uf is a Levy measure and //u^n |a;|z//(dx) < oo. Consider 
the infinitely divisible distribution A that has Levy generating triplet given 
by a = 0, Levy measure vj and 7 = — J^^iy xvf(dx) (for this notation, see 
Sato [15], pages 39 and 163); its characteristic function is 

(e izf ^ -l-izf(x))v(dx)\. 



exp 

So, by hypothesis, A = So- Hence, i>f = and thus / = 0, z^-a.e. □ 

Remark 3.6. Note that if a Levy process has finite moment of order 
k>2, then a polynomial of order [k/2] without independent term satisfies 
the conditions on the function / of the lemma. 

3.3. Kailath-Segall polynomials and finitely supported Levy measures. 

Theorem 3.7. There exists a number k>l and a family of polynomials 
{Qn(x±, ■ ■ ■ ,Xk,t),n > k}, Q n of degree n, such that 

P n {x[ 1] , . . . , X[ n) ) = Q n (X^ ,...,X { t k \t), n>k, 
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being this k the minimum number that satisfies that condition, if and only if 

l{o-^o} + #Supp(i/) = k, 
where Supp(i^) is the support of the Levy measure v. 

Proof. We first prove that the condition is sufficient. 
Case 1. Let a = and Supp(f ) = {a±, . . . ,a k }. Consider the polynomial 
of degree k + 1 , 

k k 

R(x) = x Y[(x - aj) = Cjx j + x k+1 , 

3=1 3=1 

that satisfies R(x) = 0, v-a.e. 

Denote by Lr the polynomial of order 1 in x\, . . . , x k+ \ defined by the 
coefficients of R, 

k 

L R (xi,...,x k+ i) = ^2cjXj +x k+ i, 

3=1 

and let 

c k+1 = -E[L R (x[ 1 \...,x[ h+1) )]. 
Then, the characteristic function of L R {x[ 1 ^ , x[ k+l ^)-\-Ck+it is [see (3.4)] 
\t [ (e lzR ^ - 1 - izR(x))v(dx) 



cxp 



which is equal to 1 because R(x) = 0, v-&.e. So, 

x^=-j:c 3 x^-c k+1 t 

3=1 

and it follows that 

Pk+i {X t X , . . . , x\ k+lS) ) = Q k+ i (xj. 1 ^ , . . . , X t (fc) , t) 

for some polynomial Q k +± of degree k + 1. Now, observe that for every 
n > k + 1 , the linear system 

aiSi + a\g 2 H h a\g k = a™ 



a-kQi + a\g 2 H h a\g k = a k 
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has a unique solution gi,...,g k since the determinant of the system is a 
Vandermonde one. Hence, 

k 

gjX 1 — x n = 0, u-a.e. 

i=i 

Therefore, 

where d n is defined is a similar way as before and thus 

P n {X[ l) ,...,X^ n) ) = Q n {x[ l) X[ k) ,t) Mn > k. 

Case 2. When a ^ and Supp(z^) = {a\, . . . , afe_i}, consider the polyno- 
mial of order k + 1 without independent and linear terms, 

k-l 

R(x) = x 2 Y[ (x — cij). 
j=i 

Working as in case 1, we have 

k 

Y {k+1) _ n Y { 3) n f 

X t — - ^ CjX t - C k+ it. 

Necessity of the condition. Assume that 

Pk+i{Xt l \ . . . , A"j fc+1 ^) = Q k+ i{x\ l \ . . . ,X t k ,t). 

In the left-hand side, by (1.3), the process appears to be simply just 

multiplied by l/{k + 1), so 

xf +1 Wol(xf\...,xf\t), 

where "Pol" means a polynomial in the specified variables. By Proposi- 
tion 3.2, 

k 

"(fc+i) _ V- Y U) 



xr i) =j:^ x t+ c ^ t 



and there is no linear relationship between any of the X^ , . . . , X^ . Taking 
expectations, 

k 

\ /TED / 

i=3' 



k 

-c 2 (a 2 + [ x 2 v(dx)\ - Y [ x j v{dx)+ [ x k+l v(dx). 
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Then, the characteristic function of Y^j=i c j^t^ — X[ k+l ^ +Ck+it is 1. Specif- 
ically, 

J 12.22 
exp< —-^c x tz a 



+ 1 L ( e " (Ejfc=1 ° jxj ~ xk+1) -l-iz (j2 c i xj ~ xk+1 ) ) 



If a > 0, then c\ = and, by Lemma 3.5, 



k 

x k+1 - ^2 c j xj = °> u-a.e. 

3=2 



If a = 0, then 

k 



X k+l - ^ c j xj = °> V-&.B. 

3=1 

For the second case (the first one is very similar), if the polynomial x k+1 — 
Y^!j=i c j x ^ nas onr y r <k real, nonzero, distinct roots, then Supp(^) = {a\, . . . , 
a r } and, by the sufficiency proofs above, 

(r+i) _ >^ Y (j) 
A t - 2_^CjX t +c r+ it, 

3=1 

which contradicts the assumption that there is no linear relationship between 
the X« , . . . , X( r+1 ) . So, #Supp(i/) = fc. □ 

4. Teugels polynomials associated with a Levy process. In this section, 
we will work under Nualart-Schoutens [13] conditions on the Levy mea- 
sure v, even though some definitions only need the condition that the Levy 
process has finite moments of all orders, and that can be weakened to use 
only finite moments up to a convenient order. The Nualart-Schoutens [13] 
conditions can be expressed as the existence of A > such that 

(4.1) / e x ^v(dx) <oo. 

This implies that 

/ |x|^(<ix) < oo and / \x\ n v(dx) < oo Vn > 2, 

J{\x\>l} Jr 

so Xt has moments of all orders, and the characteristic function of Xt is 
analytic. 
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Consider the measures 

^y(dx) = x 2 v{dx) and, if a > 0, j a (dx) = a 2 5o(dx) + x 2 u(dx). 

Since v has moments of all orders > 2, it follows that 7 and 7 " are finite 
measures with finite moments of all orders, and the probabilities 7/7 (R) and 
7°"/7 cr (R) have characteristic functions that are analytic in certain neighbor- 
hoods of the origin because, if we take p = A/2, then 

/ e p ^x 2 is(dx) < % ( e 2p \ x \y{dx) < 00 
7(— i,i) c P J(-i,i) c 

and also f,_ 1 ^ e p ^x 2 u{dx) < 00. This implies that the characteristic func- 
tions of 7/7 (M) and j a are determined by their moments (see Chow 
and Teicher [5], Propositions 8.4.4 and 8.4.6). 

We can construct a (finite or infinite) sequence p n (x),n > 0, of orthogonal 
monic polynomials with respect to 7 and another sequence of monic poly- 
nomials Pn(x) orthogonal with respect to 7 ". By convention, it is always the 
case that po(x) =Pq{x) = 1. 



Examples. 

1. Brownian motion. X = W . Then v = 0. The polynomials are pi(x) = x, 
p n (x) = 0,n > 2. 

2. Standard Poisson process. X = Nt — t. Then a = and v = 61, pi(x) = 
x — 1, p n {x) = 0, n > 2. 

3. Simple Levy process with 2 jump sizes, with a = 0, 

X t = ai(JVi(t) - M) + a 2 (N 2 (t) - b 2 t). 
Then v = bi5 ai + 62^2 an d 

. bia] + b 2 al 

pi{x)=x- 2 2 , 
Mi + b 2 a 2 

p 2 {x) = (x- a 1 )(x - a 2 ), 
p n (x)=0, n>3. 

4. Gamma process. (Schoutens [17]) Let {Gt,t > 0} be a Gamma process, 
that is, a Levy process such that Gt has distribution Gamma with mean t 
and scale parameter equal to 1. Consider Xt = Gt — t. Then, a = and the 
Levy measure is v(dx) = l^ x>0 j^- dx, which has infinite support. The 
sequence of orthogonal polynomials is infinite and they are the Laguerre 
polynomials Ln\x). 
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As Schoutens [17] shows, it is straightforward to obtain Pn( x ) from p n (x) 
through a family of kernels polynomials. However, in Corollary 4.3, we will 
see a useful relationship between Pn+i(%) and p n (x). 

In order to compute p n (x), write 



(4.2) 



A, 



m 2 
m 3 



m 3 
r«4 



m n+2 
m n+ 3 



m n+2 rn n+3 ■■■ m 2n +2 
where uik = J R x k ~ 2 ^/(dx) = J R x k v(dx),k > 2, and 



(4.3) 



D n {x) 



m 2 
m 3 



m 3 

ITL4 



m n+ i m n+2 
1 x 



m n+2 
m n+ 3 



m 2n +i 
x n 



There are two cases. 

1. If the support of 7 is infinite, then, for every n, A n 7^ and 

(4.4) Pn (x) = (A n _ 1 )- 1 ,D n (x) 

defines an infinite sequence of orthogonal polynomials. This follows from 
well-known facts about orthogonal polynomials; see Chihara [4], pages 51 
and 52, and Theorem 1.3.3. 

2. If v = J2k=i bj5 aj , then there are just n nonzero (7-a.e.) orthogonal poly- 
nomials po, . . . ,p n -i - The expression of the (monic) polynomial p n is very 
easily computed, as follows: 



(4.5) 



Pn(x) =(*) Y[(x-a j )=e n 1 



where 



e, 



1 Cli 



1 a r , 



1 a\ 

1 a n 
1 x 

r,n-\ 



and (*) is due to the fact that this polynomial has degree n and is identi- 
cally zero 7-a.e., so it is orthogonal to the first n orthogonal polynomials. 
For m > n, the polynomial p m (x) is also identically zero 7-a.e. and not 
unique. 
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3. When a > and v = J2k=i bj5 aj , there are also only n + 2 determinate 
polynomials p%{x), the last one being 

n 

(4.6) p n+1 (x) = x\{{x-a 3 ). 

3=1 

This expression is also deduced from the fact that this polynomial satisfies 
Pn+i = 0,7 CT -a.e. 

The orthogonal polynomials fJ^(x), or p n (x) when cr = 0, determine a 
sequence of strongly orthogonal normal martingales related to the Teugels 
martingales (see Nualart and Schoutens [13]) and that is why we call them 
the Teugels polynomials associated with X. In Section 4.3, we provide an 
explicit expression for those martingales. 

Remark 4.1. We have changed the notation of Nualart and Schoutens 
[13] and Schoutens [17] because they write p n (x) for the orthogonal polyno- 
mial of degree n — 1, and here it denotes the polynomial of degree n. 

The next theorem is a modification of the Gauss-Jacobi mechanical quadra- 
ture formula (see Szego [18], Theorems 3.4.1 and 3.4.2). 

Theorem 4.2. Assume that the Levy measure v has infinite support 
and let n > 1 be such that p n (0) ^ 0. There are then different nonzero num- 
bers a±, . . . ,a n and strictly positive numbers b\,...,b n such that the (Levy) 
measure with finite support, 

n 
k=l 

satisfies 

(4.7) / x k v(dx)= I x k v n (dx), k = 2, . . . ,2n + 1. 

Moreover, let j n (dx) = x 2 v n (dx) and ^(dx) = a 2 5o(dx) + x 2 v n (dx). Then, 
7 and 7 n (resp. 7 CT and 7^ ) have the same orthogonal polynomials up to 
order n. 

By the Gauss-Jacobi formula (Szego [18], Theorem 3.4.1), the numbers 
a±, . . . ,a n are the n different nonzero real roots of p n (x), and b\,...,b n are 
the unique solution of the compatible system 

m 2 = a\bi + + • • • + aib n 

m2n+l = «1 0\ + a 2 t>2 + ■ ■ ■ + a n ^ O n , 
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where = J* R x k u(dx), k>2. The numbers b\, ... , 6 n , called the Christoffel 
numbers, are all strictly positive (Szego [18], Theorem 3.4.2). From (4.7), it 
follows that the finite measures 7 and 7 n have the same moments of order 
0, 1, . . . , 2n — 1, denoted by 7712, ■ • ■ , rri2n+i- Hence, from the expressions (4.2), 
(4.3) and (4.4), we deduce that they have the same Teugels polynomials up 
to order n. 

Corollary 4.3. Let v he a Levy measure such that p n (0) 7^ 0. There 
are then X n and A n „i such that p^ +1 (x) = xp n (x) — A n p^(x) — A n _ip^_ 1 (x). 

Proof. The polynomials p^ +1 (x) and xp n (x) are monic and have degree 
n+ 1. The polynomial xp n (x) — p^ +1 (x) can then be written as 

11 

xp n (x)-p° +1 (x)=Y J XjP'j (x) , 
j=0 

where 

Aj = K~ l f (xp n (x) -p° n+l (x))p°(xyr(dx) 
= Kf f x Pn (x)p](x) 7 °(dx) 

and 

K d = f (pJ(z))V(<fc). 

Consider the discrete measures 7 n and 7^ of the preceding theorem. Then: 

1. the measures 7 " and 7^ have the same moments up to order 2n — 1 
and, for j = 0, . . . , n — 2, 

X j = Kr 1 xp n {x)p' T J (x)^(dx). 

J H 

2. denoting by Pj(x) [resp. (x)] the Teugels polynomials of j n (resp. 
7^), we have 

Pj(x) =pj{x) and pj(x) = pj(x), j = l,...,n, 

and [(4.6)] 

p£+l(aO = xp n {x) = xp n (x) = 0, 7n- a - e -> 
so it follows that, for j = 0, . . . , n — 2, 

Aj = iC^" 1 / xp n (x)pj(x)7^(dx) 

= ET- 1 / p£ +1 (x)pj(x) 7 £(dx) = 0. n 
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4.1. An approximating sequence of simple Levy processes. An interesting 
consequence of Theorem 4.2 is that it provides a way to construct a sequence 
of simple Levy processes that converges in the Skorohod topology to X, 
satisfying the conditions of Avram [2] in order that all variations and iterated 
integrals of the sequence converge to the variations and iterated integral of 
the limit. From the separation of zeros theorem of p n (see [4]), if p n (0) = 0, 
then p n+ i(0) ^ 0. There is then a sequence mi < < • • • / oo such that 
Pm k (0) ^ 0,Vfe. Let X k = {Xk(t),t G R} be a centered Levy process with 
diffusion coefficient a and Levy measure v mk given in Theorem 4.2. That is, 
the law of Xk is 

X k {t) = Law o-W t + aj(Nj(t) - bjt). 

J"=l 

Denote by and xi the iterated integral and the variation of order n 
of Xk, respectively, and by and X^> the iterated integral and variation 
of X, respectively. 

Theorem 4.4. Let X be a Levy process that satisfies the condition (4.1) 
and such that v has infinite support. With the above notation, for every n, 

lim P[ n) = P^ and lim X^ = X^ 

k h k K 

(both convergences in the Skorohod sense). 

Proof. By Avram [2], it suffices to prove that 

lim(Xfc, [Xk, Xk]) = (X, [X, X]) in the Skorohod sense. 

k 

Since all of the process involved are Levy process, by Jacod and Shiryaev 
[6], Corollary VII. 3. 6, it is sufficient to prove that 

\im(X k (l), [X k ,X k ]i) = (A(l), [X,X] t ) in distribution, 

and by the Cramer- Wold device, this is equivalent to proving that for every 
u,veR, 

lim(uX k (l) + v[X k ,X k ]i) = uX(l) + v[X,X]i in distribution. 

k 

From (3.4), the characteristic function of uX(l) + v[X,X]i is 
ip(z) = expj — ^u 2 z 2 a 2 + izv ^a 2 + J x 2 u(dx) 

+ t [ ( e iziux+vx2) - 1 - iz(ux + vx 2 ))u{dx) 
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From the fact that the characteristic function of X\ is analytic, it follows 
that ip(z) also is. So, it suffices to show that all cumulants of uXk(l) + 
v[Xk,Xk]i converge to the corresponding cumulants of uX(l) +v[X, X]\ 
and this is clear from the construction of v m , . □ 



4.2. The relationship between Kailath-S eg all polynomials and Teugels poly- 
nomials. 



4.2.1. Preliminary results. This subsection is purely algebraic; later, we 
will give a probabilistic interpretation of the results. First, it is convenient 
to introduce a new notation. Given a polynomial of order n, 

P(x) = c + c\x H h c n x n , 

we denote by C(P)(xi, . . . , x n +i) the polynomial of degree 1 in x\, . . . , x n+ i 
associated with the coefficients of P: 

(4.8) C(P)(xx,. . .,x n+1 ) = c xi + • • • + 

Of course, we can recover P(x) from C(P)(x±, . . . ,x n+ i): 

P(x)=£(P)(l,x,...,x n ). 

Second, we need to consider some finite-dimensional vector spaces. Let a = 
(oi, . . . , a n ), where oi, . . . ,a n are different nonzero numbers. Write 

f 

{ j=l 

n n 

X2 = y ^2,a)y j ,. . . ,x n+1 = ^2a^ +1 yj, 

for some (yi , . . . , y n ) 6 R n | . 
Sn+i 1S subspace of dimension n of R n+1 , and there is the projection 

Mn+l ra 
*■ <->n+l' 

{x\ , • • • , ► (xi , . . . , x n , u n _|_i) , 

where u n +i is computed as follows. By the Vandermonde determinant prop- 
erty, we can find (y%, . . . , y n ) £ M n such that 

n n n 

(4.9) xi = ^2ajyj, x 2 = ^ afyj, . . . , x n = ^ aJ-yj. 

j=i j=i j=i 
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We then write 



n 

a j Vr 

3=1 



Lemma 4.5. With the above notation, 

u n+ i = -C(P)(xi, . . .,x n ,0), 

where 

n 

P(x) = H(x- aj ). 

3=1 

Proof. From the expression of P(x) given in (4.5). 

1 a\ 

C(P)(x 1 ,...,x n+1 ) = 6~ 1 



1 a n 

X\ X 2 



Xn+l 



Hence, 



c(p) (n, . . . , Xn , u n+l ) = c{p) J2 a m ■ ■ • > E a l +1 yj = °- 

\j=i j=i / 

The polynomial P(x) is monic, so 

u n+ i = -£(P)(xi,...,cc n ,0). □ 

Define the polynomial of degree n + 1 , 

n 1 

J* +1 (a;i,...,z n ) =P n+ i(x 1 ,...,a: n ,it ri+ i) = E II /^jV^'K'' 

3=1 J ' 

where the summation is over all nonnegative integers hx,...,h n such that 
Yjj=i hj = n+l, [x] n is the falling factorial and y\, . . . ,y n are given in (4.9). 
This strange expression is the multiple convolution of Charlier polynomials 
aC.(yj,0),j = 1, . . . ,n. Note that when working with the polynomials, the 
variable t does not play the role of time and can be used freely according to 
our needs. 

Proposition 4.6. For every (xi,...,x n+ i) eE n+1 , 

(-1)" 



^71+1 (#1) • • • ) x n+l) ~ J n +l[%l, • • ■ , 



n+l 



£(P)(xi,...,x n+1 ), 
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where 



Equivalently, 



P(x) = U(x- aj ). 

3=1 



< — 1 ) n 

P n+1 (l, x,...,x n )- J n a +1 (l, . . . , x"- 1 ) = { —+-P(x). 

n + 1 

Proof. Simply note that P n +\ is linear in x n+ i, with coefficient (— l) n /(n + 
1) [see (1-3)], and apply Lemma 4.5. □ 

Note that this proposition is true if we replace P n +\ by another polyno- 
mial linear in the variable x n +i ; but we will see that with P n +i, it has an 
interesting probabilistic interpretation. 

In the same way, take a > and write 

f 

S n+2 = \ (xi,...,x n+2 ) £l n+2 :xi = ay + ^ j a j y j , 

n n 

X2 = ^T / Ojy j ,..., x n+2 = a l +l yj ' 
i=i i=± 



for some (y , Vi, ■ ■ ■ , Vn) G 



Consider the projection 



p(j,a 
<->n+2> 



R n+2 

(xi,...,x n+2 ) -> (xi,...,x n+ i,< +2 ), 



where 



yi,...,y n being the solution of x 2 = Ej=i tfyj, x n +i = Ej=i dj + Vj- With 
the same proof as Lemma 4.5, we have that C(P)(x 2 , ■ ■ ■ , x n+ \, u° +2 ) = 0, 
where P(x) = Yl™ =1 (x — aj). Also, note that 

C(xP)(x 1 , . . .,x n+2 ) = jC(P)(x 2 , • • -,x n+2 ). 

Define the polynomial 

J n+2 (xi,..., x n+ i) = P n+2 (xi,..., x n+ i , < +2 ) , 

which has an expression similar to J^+i with the addition of a Hermite 
polynomial H.(yQ,0). We then have the following. 
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Proposition 4.7. For every (xi,...,x n+2 ) eK n+2 , 

(-l) n+1 

P n+2 (x 1 ,...,x n+2 ) - J°' +2 (x 1 ,...,x n+1 ) = — — £(xP)(xi,...,x n+2 ), 

where 

n 

P(x) = H(x-a j ). 
3=1 

Equivalently, 

P n+2 (l, x,..., x n+1 ) - J^ a 2 (l, x, . . . , x n ) = 1 J +2 xP(x). 

4.2.2. Teugels polynomials. The propositions of the previous subsection 
can be transferred when we have a Levy measure v and the corresponding 
Teugels polynomials p n {x) and p^(x). We use Corollary 4.3 to identify these 
polynomials. 



Corollary 4.8. Fix n > 1 such that p n (0) ^ and let a±, . . . ,a n be the 
roots ofp n (x). Then, 



P n+ i(xt, . . . ,x n+1 ) - J% +1 (xi,.. -,x n )= ; - C(p n )(xi, . . .,x n+ i) 



-l) n 
n + 1 

and 



P n+2 {xi,...,x n+2 ) - J°'+ 2 (xi,...,x n+1 ) 

(-l)n+l 



2 -{^{p a n + l){xi,...,X n+2 ) 

+ K£(Pn)( x li ■ ■ -> x n+l) + \i-\£(Pn-l)( x li ■ ■ -i x n))- 

Equivalently, 

l ( — 1)™ 

fn+l(lj X,.. . , X n ) — J* +1 (l, X, . . . , X n ) = — p n (x) 

n + 1 

and 

P n+2 (l,x,...,x n+1 )-J^ 2 (l,x,...,x n ) 

_ (-l) n+1 



n + 2 



-(K+lO) + AnPn(ac) + An-lPn-l(x)). 
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4.3. Main result. The Teugels martingales {Y^ n \n > 1} (see Nualart 
and Schoutens [13]) are defined by Y t (n) = X ( t n) - E[X t (n) ], n > 1. Specifically, 

and 

y(n) = x (n) _ t f x n v ^ ^ n > 3 



(Nualart and Schoutens [13] write x[ 2) =J2 s < t (AX s ) 2 , instead of [X,X], 
as we have done; however, both definitions give the same y( 2 ).) By an or- 
thogonalization procedure, they obtain a family {H^ n \n > 1} of normal 
martingales, pairwise strongly orthogonal, that, under the hypothesis (4.1), 
generate all of L 2 (Q) by sums of iterated integrals. In order to strongly 
orthogonalize {Y^ n \n > 1}, if a > 0, they show that you can look for the 
sequence of orthogonal polynomials Pn( x ) an d take 

H^=C{pl){Y^\...,Y^), 
and the same expression with p n replacing p a n if a = 0. 

Theorem 4.9. Let X be a centered Levy process with moments of all 
orders and fix n > 1 such that the Teugels polynomial of order n, p n {x), does 
not have a zero root. Let a\,...,a n be the roots of p n {x). If a = 0, then 



(4.10) 



where 



C n = x(p n (x) -p n (0))v(dx) 

JR 

and J* +1 (X^\. . . ,x[ n) ) + ^C n t is orthogonal to Pi(X, (1) ), . . .,P n (x{ 1] , 

...,xt ] ). 

If a > 0, then 

P n+2 (X?\...,X^) = t^ {H (n + 2) +XnH (n + l )+Xn _ iH ( n)) 

n + 2 

f i ^n+l 

+ J n+2l A t J + — To ^n+lt, 

where X n and A n _i are given in Corollary 4.3 and 

Z? n+ i = a 2 p n (0) + / x 2 p n (x)v(dx). 

JR 
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Moreover, J^ 2 {x[ l \ . . . ,x[ n+1 ^) + D n+ it is orthogonal to Pi(xjp), 

Proof. Since the proof for a > is very similar to the case a = 0, we 
consider only the latter one. Formula 4.10 follows from Corollary 4.8 and 

^f +1) = £(p n )(y t (1) ,...,y/ n+1) ) 

= C( Pn ){X ( t 1 \...,X i t n+1) )-t [ x(p n (x)- Pn (0))v(dx). 



To prove the orthogonality between J^ +1 (X^\ . . . ,xj: n ') + C n t and 

/^•(A^, . . . , -2Q ) for j = 1,. . . ,n, observe that, by definition of J^+i and 
Lemma 4.5, 

ra /y(l) y(")^ , ("l)" r + _ p /y(l) y(n) T/ (n+l)\ 

J n+ll A t J + ~ + y ^nt — -Tn+1 {^i , ■ • ■ , A t , V t j, 

where 

K ( " +1) = -£(p n )(xf >,..., JfW.Oj + Ct. 

The idea of the proof is to construct a simple Levy process Z% such that, for 
r < n, 

E[p r (x«, . . . ,x«)p n+1 (x«, . . „if»y<" +1) )] 

(4.11) 

= E[P r (Z«,...,zM)P n+1 (zW,...,zf +1 ))] 

and by the orthogonality of the iterated integrals of different order, the 
expectation on the right is zero. 

With this objective, consider the Levy process (A" t , . . . , x[ n \ V^ n+1 ^) 
that has characteristic function [see (3.3)] 

<px(z) = exp(t / (jCL£t*i*'-*» + i*Pn{*)) _ i _ ix{zi _ Zn+lPn ( x) )) v (dx) 



where z = (z±, . . . ,z n+ i). 

On the other hand, let b\, . . . ,b n the Christoffel numbers corresponding to 
v given in Theorem 4.2. Define (on the same probability space or another) 

n 

Z t = Y / a j (N n (t)-b j t), 

3=1 

where N\, . . . , N n are independent Poisson processes with respective intensi- 
ties b%,...,b n . The characteristic function of (zj: , . . . , Z^ n+1 ^) is [see (3.3)] 



(z) = expji / (e 1 ^^ 1 ZjX ^ — 1 — ixz\)v n {dx" 
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where v n = Y^j=i bj$ aj - By Theorem 4.2, v and u n have the same moments 
up to order 2n + 1, both v and v n have the same first n Teugels polyno- 
mials and p n (x) =0,i/ n -a.e. Let j±, . . . , j n +x be nonnegative integers such 
that Y^uLi kjk < 2n + 1. Therefore, the (joint) cumulant of order j\ in the 
first component, order j 2 in the second component and so on, of the vectors 
{x[ l \ . . . , x[ n \ v} n+1 ^) and (Z^\ . . . , z[ n+ ) are the same. So, a polyno- 
mial up to degree 2n + 1 of (X\ l) , . . . , X[ n) , V t {n+1) ) and the same polynomial 
of (Z t (1) ,...,Z, (n+1) ) have the same expectation. In particular, for r < n, we 
have the identity (4.11). □ 

Example. A very simple example will help to interpret Theorem 4.9. 
Consider a Levy process X with a = and Levy measure v, and let a be 
the root of its Teugels polynomial of order 1, pi(x). Assume a/0 and let b 
be the solution of 

x 2 u(dx) = ba 2 . 

Let Zt = a(Nt — bt), where Nt is a Poisson process of intensity b. Zt is then 
a simple Levy process that has Levy measure v\ = b5 a . By Gauss-Jacobi 
Theorem 4.2, v and v\ have the same moments of order 2 and 3. So, Xt 
and Zt have the same cumulants of order 2 and 3, and, since both are 
centered, they have the same moments of those orders. Then, on one hand, 
,mx[ 1) )=P 2 (X u aX t ), so 

J%{X t )-\t f x 2 u{dx) = P 2 (x t ,aX t + t [ x 2 v(dx)\. 

On the other hand, 

P 2 {Z ( t 1 \ Z\ 2) ) = P 2 (Z t ,aZ t + ba 2 t). 



We then have 
E 



Pi{X t )(j%{X t ) - \t f^Mdx^j 



:E[P 1 (4 1) )P 2 (4 1) ,4 2) )] 



because Pi(X t )(J 2 (X t ) — f R x 2 v(dx)) is a product of a polynomial of de- 
gree 1 and a polynomial of order 2, in Xt, which is centered, so the expec- 
tation of that product depends only on the moments of order 2 and 3. So, 
Theorem 4.9 says that we have a decomposition 

P 2 {X[ 1) ,X[ 2) ) = J$(X t ) -\tj x 2 v(dx) - \H[ 2 \ 

JR 

such that: 

1. J 2 {X t ) - \t J R x 2 u(dx) is orthogonal to Pi(A" 4 (1) ); 
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2. J%(X t ) - \t^x 2 u{dx) = P 2 (X} 1) ,V t {2 >), where {X\ l \ v} 2) ) is a Levy pro- 
cess that has the same moments (up to order 3) as the variations {zj: , Z^) 
of the simple Levy process Z t . 

Acknowledgments. The authors wish to express their thanks to Jose An- 
tonio Carrillo and Xavier Mora for their fruitful comments about the eikonal 
equation. 
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